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Using the axially-symmetric time-dependent mean-field Gross-Pitaevskii equation we 
study the Josephson oscillation in a repulsive Bose-Einstein condensate trapped by a 
harmonic plus an one-dimensional optical-lattice potential to describe the experiments by 



O ; Cataliotti et. al. [Science 293 (2001) 843, New J. Phys. 5 (2003) 71.1]. After a study of 

^ ■ the formation of matter- wave interference upon releasing the condensate from the optical 
trap, we directly investigate the alternating atomic superfiuid Josephson current upon 

^ ! displacing the harmonic trap along the optical axis. The Josephson current is found to 

\ be disrupted upon displacing the harmonic trap through a distance greater than a critical 

^ \ distance signaling a superfiuid to a classical insulator transition in the condensate. 
OO ' 

^ ■ 1. INTRODUCTION 

o ■ 



The experimental loading of a cigar-shaped Bose-Einstein condenate (BEG) in both 
one- |l|2|3|4j and three-dimensional periodic optical-lattice potentials has allowed to 
study the quantum phase effects on a macroscopic scale such as interference of matter 
\ waves. There have been several theoretical studies on a BEC in a one- |6|7|8j and three- 
dimensional P optical-lattice potentials. The phase coherence between different sites of 
O ■ a trapped BEC on an optical lattice has been established in recent experiments |H2|3|4|5] 

through the formation of distinct interference pattern when the traps are removed. 
^ \ In the experiment on matter-wave interference from two pieces of coherent BEC, an 
5j \ interference pattern comprised of a large number of dark and bright patches is formed 
[TUj . This is similar to the well-known double-slit interference pattern in optics. As the 
number of slits in the experiment on interference of light is increased, the number of 
bright patches in the interference pattern is reduced and one has only a few prominent 
bright patches for interference of light from an optical grating with very large number of 
narrow slits. A similar phenomenon has emerged in matter-wave interference from a BEC 
in an optical-lattice trap created by a standing-wave laser field, which can be considered 
to be a large number of coherent sources of matter wave. When a BEC trapped in an 
axially-symmetric harmonic and an one-dimensional optical-lattice trap is released from 
the joint traps a definite interference pattern composed of three peaks is usually formed. 
With the increase of lattice spacing the interference pattern evolves to (2A/' + 1) peaks 
with M peaks symmetrically located in a straight line on opposite sides of a central peak. 
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Cataliotti et al. j2] prepared a BEC on a joint harmonic plus an optical-lattice trap. 
Upon displacing the harmonic trap along the optical lattice, the BEC was found to execute 
the Josephson oscillation by quantum tunneling through the optical-lattice barriers. In 
a later experiment they found that for a larger displacement of the harmonic trap 
the Josephson oscillation is disrupted. Their measurement of the Josephson oscillation 
was rather indirect and based on the existence of phase coherence and the formation of 
interference pattern upon release from the traps. They adopted the indirect procedure as 
the expansion of the BEC upon release from the traps facilitates the observation. However, 
this is not quite necessary in numerical simulation, where one can directly identify the 
Josephson oscillation and its disruption without resorting to an expansion. A preliminary 
theoretical study in this topic was based on an expansion of the BEC as in the experiment 
[7j. Here, we investigate the Josephson oscillation and its disruption based on a direct 
study of the BEC. We compare these results with experiments |2|3j as well as with previous 
results [Tj based on an expansion of the BEC upon release from the traps. 



2. MEAN-FIELD MODEL AND RESULTS 

The time-dependent BEC wave function \I'(r;r) at position r and time r is described 
by the following mean-field nonlinear Gross-Pitaevskii (GP) equation [7|8j 
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where m is the mass and the number of atoms, Q = Anh'^a/'m the strength of interaction, 
with a the scattering length. In the presence of the combined traps V^(r) = |mc(j^(p^ + 
z/^y^) + V^pt where u is the angular frequency of the harmonic trap in the radial direction 
p, uu that in the axial direction y, and V^pt is the optical-lattice potential. The axially- 
symmetric wave function can be written as \E'(r,r) = ilj{p,y,T), where < p < cxd and 
—oo<y<oo. Transforming to dimensionless variables p = \/2p/l, y = \/2y/l, t = tlu, 

I = ^h/ {muj), and ip{p, y; t) = p\/P/ V8ip{p, y; r), Eq. ((T)) becomes |7l8j 
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where n = Na/l. In terms of probability P{y,t) = 2tt JI^ dp\ip{p,y,t)\'^ / p, the normal- 
ization of the wave function is dyP{y,t) = 1. The probability P{y,t) is useful in the 
study of the formation and evolution of the interference pattern and Josephson oscillation. 

In the experiments of Cataliotti et al. |2|3j with repulsive ^^Rb atoms, the axial and 
radial trap frequencies were z/cj = 27r x 9 Hz and u; = 27r x 92 Hz, respectively. The optical 
potential created with the standing-wave laser field of wave length A = 795 nm is given by 



K)pt = VoErcos^{klz), with Er 
In terms of the dimensionless laser wave length Aq 



27r/A and Vq 
= V2A// ~ 1, 



(< 12) the strength 
ER/{hu) = 47rVA2 
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Figure 1. P{y,t) vs. y and t for the trapped BEC after the removal of combined optical 
and harmonic traps at t = for (a) c = 1, (b) c = 2.4, (c) c = 3.5 and (d) c = 4.6. The 
interference pattern has led to 3, 5, 7 and 9 peaks in these cases. 



Hence V^pt of Eq. (j21) becomes 

(3) 
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where the parameter c controls the spacing cAo/2 between the optical- lattice sites. The 
experimental condition of Cataliotti et al. |2l3j is obtained by taking c = 1. 

The GP equation Q is solved by the Crank-Nicolson method [11 . An interference 
pattern is formed by suddenly removing the combined traps at time t = on the ground- 
state solution. To study the time evolution of the system we plot in Figs. 1 (a), (b), (c) 
and (d) P{y,t) vs. y and t for c = 1,2.4,3.5 and 4.6, respectively. The variation of c of 
Eq. (jni) corresponds to a variation of the spacing between successive sites. The increase in 
c simulates an increase in the distance between the sites and a decrease in the number of 
occupied sites. In these plots one can clearly see the central condensate and the moving 
interference peak(s). As the number of occupied optical-lattice sites decreases with the 
increase of c, the interference pattern develops more and more peaks. 

We have illustrated in Fig. 1 the formation of interference pattern upon releasing the 
BEC from the joint optical and harmonic traps. The formation of interference pattern 
implies phase coherence and superfluidity in the condensate. The Josephson oscillation is 
a direct manifestation of superfluidity while the condensed atoms freely tunnel through 
the high optical lattice barriers. The absence of interference pattern after displacing 
the harmonic trap implies the loss of superfluidity and a disruption of the Josephson 
oscillation. This phenomenon represents a superfluid to a classical insulator transition and 
can be studied by a mean-field approach. The superfiuid to a quantum Mott insulator 
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Figure 2. Contour plot of P{y,t) vs. y and 
t for Vq = 2Efi after a harmonic trap dis- 
placement of 30 /im showing the Josephson 
oscillation with frequency 8.45 Hz. 
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Figure 3. Contour plot of P{y,t) vs. y 
and t for Vq = after a harmonic trap 
displacement of 70 /im demonstrating the 
disruption of the Josephson oscillation. 



transition as in Ref. can only be understood by a field-theoretic approach beyond 
mean-field theory. We study the Josephson oscillation using both approaches, e.g., upon 
releasing the BEC from the joint traps as in the previous study [7j and by following the 
condensate directly after displacing the harmonic trap. Although a free expansion is the 
only way to observe the Josephson oscillation experimentally, for numerical purpose the 
direct approach seems to be more precise and involves less computer memory and time. 

In Fig. 2 we exhibit a contour plot of P{y,t) vs. y and t for Vq = 2Eji after displacing 
the harmonic trap by a distance 30 /xm as in the experiment [2]. The Josephson sinusoidal 
oscillation around the displaced trap center at y = —30 fim is clearly visible in this plot 
from which the frequency of oscillation can be obtained reasonably accurately. However, 
when the displacement of the harmonic trap is increased beyond a critical value, the 
oscillatory motion is disrupted as shown in Fig. 3 for Vq = 5Efi and a harmonic trap 
displacement of 70 fim in agreement with experiment 3^. In this case, unlike in Fig. 
2, the condensate does not cross the center of the displaced trap at y = —70 fim. We 
performed a direct study of the Josephson oscillation for different Vq for a displacement 
of the harmonic trap below the critical value for the disruption of superfluidity. In Fig. 4 
we plot the Josephson frequencies calculated from this study as well as those calculated 
by allowing an expansion of the BEC as in Ref. [7]. In this figure we also plot ^ the 
experimental results as well as those obtained by using the tight-binding approximation. 
The three-dimensional results obtained after expansion in Ref. [7] and obtained directly 
here are in agreement with each other as well as with experiment The present direct 
results fit a smooth line and hence seem to be more accurate than the results of Ref. 
[7j. The results for tight-binding approximation [2] are slightly different from the full 
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Figure 4. Josephson frequency vs. Vq: solid 
circle with error bar - experiment f7; solid 
triangle - tight-binding approximation jSj; 
square - indirect result after expansion [7j; 
star with full line - present direct result. 



Figure 5. Axial width (Ry) of central peak 
by peak separation d vs. time spent in dis- 
placed trap: experiment — star - dis- 
placement of 30 /xm, triangle - displacement 
of 120 /im; present theory — full lines. 



three-dimensional results. 

Finally, in Fig. 5 we plot the axial width Ry of the central peak normalized to peak 
separation d for Vq = S-E/j vs. time spent in the displaced trap for displacements of 30 
Hm and 120 /xm and compare with experiment For the displacement of 30 /xm, Ry/d 
remains constant, whereas, for 120 /im, Ry/d increases to unity with time. However, the 
theoretical increase is much faster than in experiment. 

To conclude, we have provided an account of matter- wave interference, Josephson oscil- 
lation and its disruption using the three-dimensional mean-field GP equation. The results 
are in agreement with recent experiments by Cataliotti et al. j2|3j . The present results 
for Josephson frequencies are slightly different from those of tight-binding approximation. 
Further studies in three dimension are needed to understand this difference. 
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